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Abstract 

A parabolic resonance is formed when an integrable two-degrees- 
of-freedom Hamiltonian system possessing a circle of parabolic fixed 
points is perturbed. It is proved that its occurrence is generic for one 
parameter families (co-dimension one phenomenon) of near-integrable, 
two d.o.f. Hamiltonian systems. Numerical experiments indicate that 
the motion near a parabolic resonance exhibits a new type of chaotic 
behavior which includes instabilities in some directions and long trap- 
ping times in others. Moreover, in a degenerate case, near a flat 
parabolic resonance, large scale instabilities appear. A model aris- 
ing from an atmospherical study is shown to exhibit flat parabolic 
resonance. This supplies a simple mechanism for the transport of par- 
ticles with small (i.e. atmospherically relevant) initial velocities from 
the vicinity of the equator to high latitudes. A modification of the 
model which allows the development of atmospherical jets unfolds the 
degeneracy, yet traces of the fiat instabilities are clearly observed. 
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1 Introduction 

The motion of n particles under conservative forces (e.g. the plan- 
ets in the solar system, coupled frictionless pendulums, point vor- 
tices etc.) is described by an n-degrees-of-freedom Hamiltonian 
system. This motion may be ordered or chaotic, depending on 
the Hamiltonian H{x,v) (the energy), the dimension of the system 
and the initial position and velocities of the particles. Consider 
a Hamiltonian flow for which all solutions are ordered. A funda- 
mental question is what will be the solution's fate if the system is 
changed just by a small amount (e.g. the efl"ect of a meteor on the 
solar system). In particular, consider solutions which are bounded 
and are "roughly stable" in the unperturbed system (there exists 
a flnite size neighborhood in phase space which stays "near" these 
solutions forever). Can solutions in this neighborhood become un- 
bounded or considerably difl"erent if the system is changed by lit- 
tle? In this paper we discover a new low dimensional mechanism 
(n = 2) creating such a dramatic efl"ect on solutions. We prove it 
is a mechanism which is physically common, and demonstrate its 
consequences on a simple model for the motion of weather balloons 
in the atmosphere. 

The appearance of homoclinic loops in integrable Hamiltonian systems 
provides a mechanism by which by small perturbation integrability is de- 
stroyed and new complex dynamics arises. This phenomenon has been prin- 
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cipally investigated in the context of one-and-a-half d.o.f. Hamiltonian sys- 
tem. Yet, even with regards to this simplest setting for which homoclinic 
chaos exists, important open questions linger, see for example papers in this 
issue by Klafter, Meiss, Melnikov and Zaslavsky and references therein. 

What is the structure of near-integrable two d.o.f. Hamiltonian system? 
Generically, most phase space is foliated by two-dimensional tori which divide 
the energy surfaces. Substantially different dynamics appear near singular 
points or curves on the energy surface - either fixed points or invariant circles 
(or lines in the non-compact case) of the fiow [15]. When isolated unstable 
fixed points exist inherently three dimensional phenomena arise, and these 
are not yet fully understood. See Lerman and Umanskii [15]for the classifica- 
tion of integrable systems and Koltsova and Lerman [16]and Grotta-Ragazzo 
[14]for the analysis of chaotic system with a homoclinic loop to a saddle- 
center. 

Consider near-integrable Hamiltonians of the form: 

H{x, V, D, 9; e) = Ho{x, v, D)+eHi{x, v, D, 9; e), (x, v, D) G i?^ ^ G T, 

(1-0.1) 

Such a form arises from a general Hamiltonian H{x, j/, w] e) of a near inte- 
grable fiow by using one of the constants of motion of the unperturbed fiow 
(e = 0), D[x^v^y^w) as a variable, and defining 9 as its conjugate variable 
(e.g. D is the angular momentum). Consider a domain U for which the 
transformation (x^v^y^w) — )• [x^v^D^9) is non-singular^. In particular, in 
[/, a fixed point in the (x, v) plane, qj = (xj, uj, Dj): 

V^,,,)Hoixf,Vf,Df) = (1.0.2) 

corresponds to an invariant circle in the four-dimensional phase space (and 
not to an isolated fixed point). This circle may be composed of a fam- 
ily of periodic orbits (when 9[qf) ^ 0) or of a circle of fixed points (when 
9{xf^ uj, Df) = 0). An important lesson from Lerman and Umanskii [15]work 
is that phenomena which are co-dimension one in near-integrable one-and- 
a-half d.o.f. Hamiltonian systems may become generic in integrable two 
d.o.f. Hamiltonian systems; In the (x^v) plane, qj may generically be hy- 
perbolic or elliptic, and generically, there exists Dj = Dp values for which 

"'^i.e. where the surfaces spanned by the Poisson action of the integrable system [15] are 
at least one-dimensional 
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it is parabolic [15]. Similarily, generically, there exists Dj = Dr values lor 
which 6{x f ,v f , D f) = 0. These simple observations are the basis lor the 
prool that parabolic circles ol fixed points (i.e. Dj = Dr = Dp) are generic 
in one-parameter lamily ol integrable two d.o.l. Hamiltonian systems, see 
section 3. 

The behavior near the invariant circles under small perturbations depends 
on their stability (i.e. the stability ol qj in the (x,u) plane), on the rotation 
rate on them (i.e. 6{qf )), and the rotation rate on orbits which are homoclinic 
to them (il such exist). II qj is elliptic, stability reigns, though one may 
obtain resonant or nonresonant behavior depending on the frequency ratios. 
In particular, il qj is an elliptic circle ol fixed points, the resonance creates 
localized structures in 6. When qj is hyperbolic, and the energy surlaces are 
compact, under small perturbation several types ol homoclinic tangles may 
appear. It is conjectured that only one type, the reducible homoclinic tangle, 
is topological (orbital) conjugate to the one appearing in the one-and-a-hall 
d.o.l. Hamiltonian system case. 

The reducible homoclinic tangle appears when, in the integrable limit, 
there exists a homoclinic loop to a periodic orbit, and the motion along the 
homoclinic loop and its lat neighborhood preserves the angular direction ol 
the periodic orbit (i.e. |^| > a > along the periodic orbits and the homo- 
clinic loop). In this global Poincare map transverse to the periodic 
orbit is well defined in a neighborhood ol the homoclinic loop and the fiow 
is conjugate to a two-dimensional symplectic map with a homoclinic tangle. 
For separable two d.o.l. Hamiltonian systems homoclinic loops to periodic 
orbits ol bounded periods always satisly this requirement as in Holmes and 
Marsden [10]. In non-separable systems, it is satisfied in regions in phase 
space lor which one ol the angle variables varies monotonically in t along 
the unperturbed separatrices as in Holmes and Marsden[ll]. Even in this 
simplest setting there are some differences between two d.o.l. Hamiltonian 
systems and one-and-a-hall d.o.l. Hamiltonian systems when one considers 
behavior ol physical ensembles ol initial conditions [1, 2] or relations between, 
lor example residence time and escape rates in the Poincare section. These 
differences seems to be physically important yet mathematically - technical. 

A second type ol homoclinic tangle [lumpy homoclinic tangle) arises when 
small perturbations are applied to systems lor which non-monotonicity in 
the angle occurs along the unperturbed homoclinic solutions. Then, non- 
uniformity in the angle variable ol the chaotic zone is created [1]. Beyond 
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the differences in terms of tlie observed non-uniformity between tlie first and 
second types, it is conjectured tliat tliis case is dynamicaiiy distinct from 
tlie reducible liomoclinic tangle, at least for small yet finite perturbation 
parameter. 

A third type of homoclinic tangles, the homoclinic hyperbolic resonance , 
occurs when a system containing a normally hyperbolic circle of fixed points 
with a family of heteroclinic orbits connecting points on the circle is per- 
turbed, see Kovacic and Wiggins [8], Haller [9] and Kaper and Kovacic [7]. 
Using ideas as in Lerman and Umanskii [15], it is proved in section 3 that the 
occurrence of a circle of fixed points which is normally hyperbolic whithin 
the energy surface is generic in the class of integrable two d.o.f. Hamiltonian 
systems. Thus, if the energy surfaces are compact, the scenario of homo- 
clinic hyperbolic resonance appears generically in near-integrable two d.o.f. 
Hamiltonian systems. 

In this paper a new type of chaotic behavior is introduced - the homoclinic- 
parabolic resonance. Parabolic invariant circles appear generically in inte- 
grable non-separable two d.o.f. Hamiltonian system [15]. When there exists 
periodic motion on the invariant circles the effect of small perturbation on 
initial conditions in its neighborhood seems to be insignificant (in fact the 
separatrices splitting is exponentially small in the distance from the bifurca- 
tion point [12]). In particular, initial conditions starting close to the invariant 
circle stay close to it - the invariant circle is "roughly stable". This changes 
dramatically when the invariant circle is a circle of fixed points in the four 
dimensional phase space, see figure 2.0.1 - figure 2.0.5. First, it is proved that 
this event is not so exotic - it is a co-dimension one phenomena for integrable 
two d.o.f. Hamiltonian system. Then, it is demonstrated numerically that 
small perturbations lead to dramatically different dynamics than observed 
in the other types of homoclinic tangles or the other types of resonances. It 
is a combination of the hyperbolic homoclinic chaos, the localization of the 
elliptic resonances and instabilities formed by sliding along the newly formed 
elliptic circles as seen in slow passage through bifurcations [13]. These ef- 
fects are even more pronounced when additional degeneracy appears - the 
fiat parabolic resonance case. Here, in the unperturbed system the parabolic 
invariant circle of fixed points is a part of a family of invariant circles of fixed 
points all of which belonging to the same energy level set Hq = h. Namely, on 
this energy surface the non-degeneracy conditions for the iso-energetic KAM 
theory fails to exist and instabilities reign. This may be an example for a 
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two-d.o.f. equivalent of the stochastic-webs [18] which appear in one-and-a- 
half d.o.f. Hamiltonian system(global instabilities will appear if non-compact 
energy surfaces are considered). 

This study has been initiated while investigating a simple model for parti- 
cles (weather balloons or floats) transport in the atmosphere or ocean. With 
Y. Dvorkin and N. Paldor [1], we have analyzed the simplest possible two- 
dimensional motion on a rotating sphere, due only to Coriolis force (i.e. in 
the absence of any body force, see Paldor and Killworth [6]), perturbed by 
the inclusion of a zonally travelling pressure wave which perturbs the geopo- 
tential surface from its mean spherical shape, see Paldor and Boss [5]. In 
[1] we developed tools for delineating the phase space of non-separable two 
d.o.f. Hamiltonian system to regions in which the different types of homo- 
clinic structures appear. Here, it is shown that this model exhibits a flat 
parabolic resonance when the perturbation is of the form of a standing pres- 
sure wave. Adding a latitude dependent pressure gradient, as was suggested 
in [1], removes this degeneracy. If the pressure gradient is comparable to 
reasonable atmospherical velocities (the atmospherical relevant case) then 
strong instabilities induced by the nearby flat parabolic resonance are ob- 
served. 

The parabolic resonance is a low co-dimension phenomenon hence it is 
expected to appear in many other applications. For example, the motion of a 
particle in a central field [4, 11] with effective potential with at least one max- 
imum is an integrable two d.o.f. Hamiltonian system, which, when perturbed 
by small amplitude travelling wave [2] (in the angle direction) may exhibit 
a parabolic resonance for some isolated wave speeds (this is a variation on 
the example of Holmes and Marsden [11] who considered this problem in the 
hyperbolic regime and with standing wave perturbation). Another example 
of similar form is the whirling pendulum perturbed by a travelling wave per- 
turbation. Notice that these examples (and the atmospherical model) are 
formally two-and-a-half d.o.f. problems (since time is explicitly introduced 
in the travelling wave perturbation), but, as is demonstrated in section 2 
and discussed in section 3, this class of systems can be easily reduced to 
autonomous two d.o.f. Hamiltonian system. Finally, the example which 
motivated the works on the homoclinic hyperbolic resonance (see Kovacic 
and Wiggins [8]), namely that of a two-mode truncation of the damped and 
driven Sine-Gordon equation (see McLaughlin paper in this issue and refer- 
ences therein) also exhibits parabolic resonance for some parameter values. 
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This paper is organized as follows: In section 2 the simple quasi-inertial 
model for Lagrangian motion in the atmosphere is recalled and studied nu- 
merically near its parabolic resonances. In section 3 it is established that the 
occurrence of a parabolic resonance is generic for near-integrable two d.o.f. 
Hamiltonian system depending on one parameter. Then several physically 
typical forms of two d.o.f. Hamiltonian system (such as separable systems) 
are considered, and their compliance with the genericity assumptions is ex- 
amined. Section 4 is devoted for discussion. 



2 Model for particles' motion in the Atmo- 
sphere 

The motion of a particle on a rotating sphere subject to a conservative travel- 
ling wave perturbation serves as a simple model for the motion of particles in 
the atmosphere [6, 5, 1]. It was recently suggested that the inclusion of lati- 
tude dependent pressure gradient may be incorporated into this model to sim- 
ulate the appearance of jets in the atmosphere [1]. Using polar co-ordinates, 
the non-dimensional Lagrangian momentum equations for the eastward and 
northward velocity components (u, v) and the rate of change of the longitude 
and latitude coordinates (A, (f>) in the presence of a zonally travelling pressure 
wave and a latitude dependent pressure gradient are given by: 

dX u 
dt cos cf) 

du -It ^ \ 7 ^(^) 

— = usm^flH -) — ke -cosikX — crt). 

dt cos (f) cos (f) 

dcf) 

dt ~ ^ 

dv u 

— = -usm(l)(l + -) - B'U) - eA'U)sm(kX - at) 

dt cos (f) 

(^G [-^,^], AG [0,27r], u.veR'. 

B[(f>) and A[(f>) represent the latitude dependent amplitudes of the pressure 
gradient and the pressure wave respectively. For simplicity, both are assumed 
to be even in (p. For e = 0, the system is autonomous, A is uncoupled and 



(2.0.3) 
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(2.0.3) has two integrals of motion, corresponding to (twice) the angular 
momentum D and the energy E: 

D = cos (^(cos + 2m) (2.0.4) 

E = ^{u^ + v^) + B {(!>). (2.0.5) 

In [1] the structure of (2.0.3) is studied for B[(f>) = for general values 
of the parameters. By constructing the colored energy surfaces H = h and 
the corresponding bifurcation diagrams of the energy-momentum map (in- 
dicating singular curves, their stability and the monotonicity properties of 
6[t) in the (H^D) plane), the phase space is delineated to regions in which 
inherently three-dimensional dynamics appears and to regions in which the 
system is reducible to one-and-a-half d.o.f. Hamiltonian system [1]. Critical 
parameter values for which the dynamics changes qualitatively are found by 
this analysis. Specifically, the case of standing waves (cr = 0) has been iden- 
tified as a degenerate case. It is shown below that for B[(f>) = 0, cr = 0, a 
fiat parabolic resonance occurs, and that by appropriate choice of B[(f>) the 
resonance may be made non-degenerate. Moreover, the behavior near elliptic 
and hyperbolic resonances is compared with that appearing near a fiat and 
generic parabolic resonances. 

Consider the case k ^ 0^ and define: 



e = ^{\-ct), c=^. (2.0.6) 



Then (2.0.3) becomes: 



de _ I D 1 1 

~di ~ 4cos2(^ ~ 2^^^ 2^ 
dD 

— = -2keA{(f)) cos{2ke) 
dcf) 

dt ~ ^ 

§ = lsM2^)il-^)-B'i^)-eA'i^)sM2k9) 
dt o cos^ (p 

which is an autonomous system with Hamiltonian 

H{<f>,v,e,D-k,c,e) = + - cos <f>f -^D + B{<f>) + eA{<f>)sin{2ke) 

2 o cos (p 2 
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k^O, < c, e > 0, sign(A;) = sign(c). (2.0.7) 

Since A[(f>) and B[(f>) are assumed to be even^, the cylinder (f> = v = is 
invariant for all e. On it the equations for (D, 6) correspond to the pendulum 
equations with 2k islands of width y^8eA(0) centered at D = D,r{c) = 1 -\- 2c, 
and rotational orbits filling the rest of the cylinder [1]. The behavior near 
this cylinder in the four-dimensional phase space depends on its stability 
properties in its normal ( [(p^v)) direction. 

Consider the unperturbed fiow in the v) plane. The origin is hyperbolic 



for |_D| < y 1 — 4i?"(0), having two symmetrical homoclinic orbits extending 
to ib(^/i,0) {(ph ~ arccos(_D) for B[(f>) <^ 1). It is parabolic at D = Dp = 



1 — 4i?"(0) and elliptic for \D\ > Dp. Parabolic resonance arises when 
Dp Dr{c), namely for c values near: 

1 . 1 



^[Dp-l) = -{,Jl-iB"{Q)-l). (2.0.8) 



For \D\ < Dp there are two elliptic fixed points (ib^e//,0) at 

Df{<peu) = COS^{(pell] 



\ sm{2<peu)' 



(2.0.9) 



with the natural frequency in the ^, v plane of nearly sin (peii (for small 
B'[(f>)^ B"[(f>))^ and with the frequency in the longitude direction of 



Notice that for B'{(f)) < 1 

It follows that when B[(f>eu) = (more generally if B'[(f>) = asin(2^)) and 
c = c„ = ( 2a) a degenerate situation occurs by which 



^Thus this assumption of symmetry is convenient, yet it seems to be inessential for 
what follows. 
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^\{4>eii,o,e,Dg{4>eii)) = = 0. Namely, lor each (p^u G (-7r/2, 7r/2), the circle 

(^e//, 0, ^, -DjJ~(^e//)) consists ol fixed points in the lour dimensional phase 
space. Moreover, lor e = all these circles ol fixed points belong to the same 
energy surlace H = 0. In section 3 it is established that this flat parabolic 
resonance is not generic. 

Numerical simulations ol (2.0.7) are performed using DSTOOL [17] with 
the Bulirsch-Stoer integrator. The Hamiltonian value is monitored so that at 
least 5 significant digits are unaltered in each computation. In all figures k = 
3,A(^) = cos"^(^), = |/3cos"^(^), and the other parameters, including 

/3 = — i?"(0), are varied as indicated in the captions. 




Figure 2.0.1: Flat parabolic resonance. 
One trajectory shown in the a) (^, 6) and b) (D, 6) planes. 
€ = 0.00025, c = 0,\\{<Po,vo,uo)\\ = O{10-'), B{<P) = 0. 

figure 2.0.1 demonstrates that strong instabilities appear near a fiat parabolic 
resonance. In this figure one trajectory starting close to the fiat parabolic res- 
onance at Dq = 1 is shown. The trajectory which is positioned initially close 
to the equator with very small initial velocities (10 cm/second = Uq = lO""*) 
travels to latitudes as high as 60° (similar behavior is observed for initial 
westward velocities ol up to Im/second). Moreover, the motion is restricted 
for a long period ol time to one cell ol the resonance band (zonally localized) 
and then jumps to a different cell. In the (^, D) plane it is seen that the mo- 
tion for D > Dp follows closely the resonant motion ol the invariant cylinder 
(f> = V = whereas the instabilities are created for D < Dp when the motion 
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descends along resonant elliptic points up to very small D values [D = 0.1 is 
seen in the figure) creating the instabilities that are seen in the 6) plane. 
As —Mo is increased (above —Uq 0.01) the descendent along these elliptic 
points disappears, and as it is increased further (above —Uq 0.05) localized 
structures disappear. Therefore, a mechanism for the transport of particles 
launched near the equator with ||(mo,^^o)|| ^ 1 to high altitudes has been 
found, a phenomena that has been observed in atmospherical experiments 
[3]. Curiously, it appears that the instabilities develop only for sufficiently 
small initial velocities. Typical initial conditions seem to be trapped in one 
zone for a long periods of time (e.g. figure 2.0.2). 



Figure 2.0.2: Near a fiat resonance. 
The same initial condition is evolved with different wave speeds: 
€ = 10-\B{<P) = 0, ((^0, vo, 00, uo) = (10-^ 10-^ 2.2, llQ-''). 
a) c = 0.1, b)c = 0.01, c) c = 0.0001. 

In figure 2.0.2 the behavior near the parabolic resonance is examined as 
the wave velocity c is varied: the same initial condition (with Dq = 1) is 
integrated three times with c = 0.1,0.01 and c = 0.0001 respectively. It is 
seen that for c values sufficiently far from the equator vicinity is stable 
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whereas when c = 0.0001 the instabihty and the locahzation of the nearby 
flat resonance is clearly seen. 




Figure 2.0.3: Nearly flat parabolic resonance. 
One trajectory shown in the a) {(p^O) and b) (D^O) planes, 
e = 0.00025, c = 0.029, B"{0) = -0.03, (Vo, ^^o) = 0(10-^), uq = 0.00042. 

If gf^24>"ii) 7^ const then the flat degeneracy is removed. The effect of 
removing the flat degeneracy is examined by varying /3 = B"[0). In figure 
2.0.3 the nearly fiat behavior is examined, with /3 = 0.03 (cp = 0.0291, Dp = 
1.058), showing that the fiat instabilities prevail. 

Figures 2.0.4-2.0.5 show the non-degenerate case with /3 = 0.3 near Dp = 
1.48 and 0.24. In figure 2.0.4 a trajectory starting very 

close to the invariant cylinder exhibits strong zonal localization and a mixed 
behavior of motion near the elliptic points, near the separatrices and near the 
separatrices of the resonance zone, figure 2.0.5 shows the behavior near the 
resonance when both c is slightly shifted from Cp and the initial condition 
of the trajectory is not extremely close to the invariant cylinder as in the 
previous figure. The chaotic trajectory is presented together with two other 
trajectories which are on the invariant cylinder (^o = t^o = 0) to illustrate 
the shift of the localized structure from its base. 

We would like to contrast the behavior observed in these figures, with 
that observed near elliptic resonance and near hyperbolic resonance; In the 
elliptic-resonance case near the resonance zone the D variable perform ex- 
cursions of Oi^^/e). However, in the plane these excursions are nearly 
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Figure 2.0.4: Parabolic resonance. 
One trajectory shown in the a) 6) and b) (D, 6) planes, 
e = 0.00025, c = 0.24, B"{0) = -0.3, ((^o, ^^o, ^o, uq) = {10~^\ 10~^\ 0, 0.24). 



unnoticeable - the motion continues to exist mostly on invariant tori. Some 
of these tori correspond to oscillatory motion in ^, and small changes in (^, v) 
usually do not alter the oscillatory/rotational behavior in ^, unless the ini- 
tial conditions are in the extreme (exponentially small?) close vicinity of the 
separatrices created by the resonance. Thus, in the elliptic-resonance case 
one observes either regular zonally localized motion or regular motion trav- 
elling around the globe and no transitions between the two type of motions 
is detected by simple numerical experiments (for small e). 

The behavior in the vicinity of a hyperbolic-resonance is mostly influenced 
by the instability in the hyperbolic direction. Thus, the 0{^/€) oscillations in 
D are coupled to the unstable chaotic motion in the (^, v) plane as is shown in 
figure 2.0.6 and figure 2.0.7. The extend in latitude to which chaotic trajec- 
tories reach is approximately the same as the latitude reached by an unper- 
turbed orbit starting with the same initial velocity ^^^x — ^max + ^(a/^) (^^^ 
figures). The perturbed motion near the invariant cylinder is chaotic w.r.t. 
north-south motion and it has a zonal structure scarred by the resonance on 
the invariant cylinder (see figure 2.0.6 and figure 2.0.7). Localized motion in 
6 exists on the invariant cylinder inside the resonance zone. Applying the 
general theory of Haller [9] for the case B[(f>) = O^k > 3, it follows that 
for c G (cfc = — sin^(^),0) homoclinic orbits to oscillatory periodic orbits 
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Figure 2.0.5: Near a parabolic resonance. 
Three trajectories are shown in the a) 6) and b) (D, 6) planes, 
e = 0.00025, c = 0.2(cp = 0.24), 5"(0) = -0.3, 
Two trajectories are on the invariant cylinder (f> = v = 
Third trajectory has (^o, ^^o, ^'o, Uq) (10~^, 10"^, 0, 0.2). 



exist [1], i.e. that localized homoclinic solutions exist (and may be detected 
numerically by devicing the correct scheme [8, 9, 7]). S -0.066, 
this range of c values is also influenced by the flat parabolic resonance, and 
so the behavior described above for the near-parabolic-resonance regime flts 
this regime as well. It is unclear yet if this is true in the general case. Taking 
smaller c values (c = —0.25), one flnds that trajectories progress along the 
longitude coordinate (with known "speed" [1, 9]), thus no zonally localized 
motion is detected outside the invariant cylinder (f> = v = 0. 

Generic parabolic resonance combines the behavior of the hyperbolic and 
elliptic regimes - one observes the localized structures due to the elliptic 
nature of the trajectory for D > Dp and then the hyperbolic nature which 
induces instabilities eventually breaking this localized structures for D < Dp. 
Moreover, another dominant mechanism for the instability is induced by 
the elliptic orbits created at D = Dp, similar to trajectories passing slowly 
through a pitchfork bifurcation [13]. Combining these three mechanisms one 
obtains a motion which has chaotic nature w.r.t. to zonal localization (long 
trapping periods and then a jump to a different island) and instabilities 
in the latitude extent of the trajectories. Initial conditions which in the 
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Figure 2.0.6: Hyperbolic resonance, = 0. 
One trajectory shown in the a) 6) and b) (D, 6) planes, 
e = 0.00025, c = -0.25, B{(P) = 0, ((^o, ^^o) = 0(10-^), uq = -0.25. 

unperturbed case remain in the vicinity of the equator (starting with D > 
Dp) may reach latitudes as high as 25° in the non-degenerate case (figure 
2.0.4) and much higher latitudes in the fiat and nearly fiat cases (figure 2.0.1 
and figure 2.0.3). Changing the parameters to values which are close to those 
producing parabolic resonances, show that much of the properties described 
above are still present - though the zonal localization occurs for lower D 
values as the distance from the invariant cylinder is increased (figure 2.0.5). 

3 Genericity of parabolic resonances. 

Consider a near-integrable two-degrees-of-freedom Hamiltonian fiow depend- 
ing on a parameter c: 

H{x,v,D,e;c) = Ho{x,v,D;c) + eHi{x,v,D,e;e), (3.0.12) 

{x,v,D) e R^, e eT, R. 

Such a form arises from a general Hamiltonian H{x,v,y,w]c,e) of a near 
integrable fiow by using one of the constants of motion of the unperturbed 
fiow, D[x^p) as a variable, and defining 6 as its conjugate variable. Consider 
a domain U for which the transformation (x^v^y^w) — )• (x^v^D^O) is non- 
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Figure 2.0.7: Hyperbolic resonance, = 1-8. 
One trajectory shown in the a) 6) and b) (D, 6) planes, 
e = 0.00025, c = -0.25, B{(P) = 0, ((^o, ^^o) = 0(10-^), uq = -0.25. 

singular"^. In particular, in [/, a fixed point in the (x^v) plane corresponds 
to either a periodic orbit or to a circle of fixed points in the four- dimensional 
phase space (and not to an isolated fixed point). 

Generally, one would expect parabolic resonances to be of at least co- 
dimension 2 in this setting (e.g. one condition for parabolicity and another 
for vanishing frequency on the invariant circle). Loosely, in view of Lerman 
and Umanskii [15]work, one parameter may be replaced by the constant of 
the unperturbed motion, D. 

It is proved below that transverse intersections of various manifolds in 
the four dimensional space (x, D, c) imply the genericity of the hyperbolic 
(resp. parabolic) resonances in near-integrable, (resp. one parameter family 
of near-integrable) two d.o.f. Hamiltonian systems. 

Surfaces of fixed points of the two-dimensional fiow in the (z, v) plane are 
determined by: 

V(^,„)i7o(9/) = 0, qj = [xf, Vf, Df, Cf), (3.0.13) 

hence these are two dimensional surfaces. Denote the union of these surface 
by F and denote points on F hy qj. Generically the fixed points qj ^ F are 

^i.e. where the surfaces spanned by the Poisson action of the integrable system [15] are 
at least one-dimensional 
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either elliptic or hyperbolic. Parabolic points on F are determined by the 
additional condition 



{d^ H^{x,v,D -c))^ 
(dx, dv) 



det ^ r I.. = o (3.0.14) 



hence correspond to a one dimensional curve P C F. 

Define the three-dimensional surface T consisting of points qx at which 
^ = 0- 

^l,T=0 (3.0.15) 

Generically, T intersects F transversely along a one dimensional curve R of 
fixed points qn R for which ^ = 0, namely q^ correspond to a circle of 
fixed point in the full phase space. Generically the fixed point on R are 
composed of intervals of elliptic fixed points and intervals of hyperbolic fixed 
points. Since generically, 

-Q^Un^O (3.0.16) 

the curve R intersect the three dimensional surface c = const transversely at 
a point (or several points), hence, we have established: 

Lemma 1: The occurrence of a normally hyperbolic (resp. elliptic) circle of 
fixed points is generic in integrable two d.o.f. Hamiltonian system. 

Corollary 2: The occurrence of a hyperbolic (elliptic) resonance is generic 
in near-integrable two d.o.f. Hamiltonian systems. 

Moreover, generically the three dimensional surface T intersects the one 
dimensional curve of parabolic fixed points P at isolated points qpR^ which 
are circles of parabolic fixed points in the four dimensional space. Thus, it 
has been established that: 

Lemma 3: The occurrence of a parabolic circle of fixed points is generic in a 
one-parameter family of integrable two d.o.f. Hamiltonian systems. 

Corollary 4: The occurrence of a parabolic resonance is generic in a one- 
parameter family of near-integrable two d.o.f. Hamiltonian systems. 
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Finally, notice that the surface of constant energy containing the parabolic 
circle of fixed points Ho[q) = Ho{qpR) defines a three dimensional surface 
HpR, which generically intersects the one dimensional curve of circles of fixed 
points, i?, at isolated points, and the two-dimensional surface F along a curve 
FpR. Hence, generically, on its energy surface, qpp is an isolated circle of 
fixed points belonging to a one dimensional family of invariant circles (Fpp) 
on which f 7^ 0. 

The tangential parabolic resonance occurs when the intersection of Hpp 
with R is not transverse at qpp^ namely Fpp and R are tangent at qpp. 
Clearly this is a co-dimension two phenomena. Then, there is an (infinites- 
imal) band of circles of fixed points on Hpp. A flat parabolic resonance 
appears when the two curves R and Fpp coincide. If along Fpp c remains 
constant, then a whole line of invariant circles of fixed points belongs to the 
same energy surface for the same parameter value - a degenerate situation 
in which no barriers to transport exist under small perturbations. The at- 
mospherical model described in section 2, exhibits, when B[(f>) = 0, this 
degenerate fiat parabolic resonance. 

The form (3.0.12) arises naturally in applications, four of which were 
mentioned in the introduction. In the applications, the Hamiltonian Ho[x^ D; c) 
may appear in some special form. The implications on the non-degeneracy 
conditions mentioned above are examined for three physically common cases: 



3.1 The separable case. 

The separable case arises naturally when two one-d.o.f.-systems are weakly 
coupled. Then, the transformation (x^v^y^w) — )• (x^v^D^O) is applied near 
invariant circle in the (y^w) plane, namely (D^O) are the action-angle co- 
ordinates of the (y^w) unperturbed fiow. Moreover, it is natural in this 
case to assume that the parameter c appears only in one of the one d.o.f. 
Hamiltonian systems. Hence we assume 

Ho{x, V, D- c) = H^{x, V- ci) + H^{D- C2), (3.1.1) 

or, more explicitly 

V,,,^ = 0, (3.1.2) 
and consider the two cases Ci = const and C2 = const separately. 
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Fixing C2, i.e. letting 

does not alter the previous results; Note that the curves corresponding to 
parabolic invariant circles, P, (which generically occur for isolated Ci values), 
are just straight lines of the form {qp = (xp, Up, Cip, D), D G [-Di, -D2]}- Since 
in this case: 

T = {qT: '^^^^§^Ur = 0} = {qr : Dt = Dn} (3.1.4) 

it follows that P intersects T transversely. Similarly, the curve R lies in T and 
generically intersects the surface c = const transversely. It follows that the 
previous Lemmas and corollaries (1-4) are valid with the words "two d.o.f. 
Hamiltonian system" replaced by "two d.o.f. Hamiltonian systems satisfying 
(3.1.2) and (3.1.3)". 



Fixing Ci, corresponds to: 



dH 

V,,,— ^ = (3.1.5) 
oc 



and produces different behavior; In this case, generically, parabolic invariant 
circles do not appear. The two dimensional surface of fixed points F, is 
parallel to the (D, c) plane, and on it there exists the curve R on which 
circles of fixed points live. Thus for this case Lemma 1 and Corollary 2 apply 
yet Lemma 3 and Corollary 4 are not valid; in fact: 

Lemma 5: A necessary condition for the generic existence of parabolic 
circle of fixed points in a one parameter family of integrable two d.o.f. Hamil- 
tonian system is 

V,,,— ^ + V,,,— ^^0. 3.1.6 
oc oD 

3.2 The travelling wave (TW) case. 

Consider unperturbed Hamiltonians of the form: 

H^{x, V, D- c) = H'^ix, V, D) + H^{D- c), (3.2.1) 
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namely, 

V,,,— ^ = 0, V,,,— ^ ^ 0. 3.2.2 
oc oD 

These are called travelling wave Hamiltonians by the following motivation; 
consider an integrable two d.o.f. Hamiltonian system with canonical coor- 
dinates (a;,u,_D,A), which is perturbed by a travelling wave perturbation in 
A: 

H{x, V, D, A) = Ho{x, V, D) + eHi{x, v, D, kX - at). (3.2.3) 

Formally this Hamiltonian corresponds to a two-and-a-half d.o.f. Hamilto- 
nian system. Using the trivial change of variables (for k ^ 0) 

e = X-yt = X-ct (3.2.4) 
k 

it is reduced to an autonomous two d.o.f. Hamiltonian system of the (TW) 
form with: 

H{x, V, D, 9) = Ho{x, V, D)-cD + eHi{x, v, D, kO). (3.2.5) 

Indeed, 

a TT 

e = A-c=— (3.2.6) 

D = _ JH,{x -^D^a) ^JH^ 
da oO 

(3.2.8) 

Three of the four examples mentioned in the introduction are of this form. 

In this case the parabolic invariant circles occur for isolated fixed D values 
independent of c, so the curves P are straight lines parallel to the c axis. The 
curves R of circles of fixed points are as in the general case, hence Lemma 
1-4 are valid for the TW case as well. 



3.3 The natural mechanical case. 

If Hq is constructed from a family of one-degree-of-freedom systems with an 
effective potential which depends on D: 

Ho{x,v,D;c) = ^v^ + V{x,D;c) (3.3.1) 
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then qpR = [xpR^O^ DpR^ cpr)^ and the conditions (3.0.13), (3.0.14), (3.0.15) 
become: 



dx '"^^ 



d 



2 \<1PR 



(3.3.2) 



Generically this form obeys the assumed transversahty of the manifolds as 
in the general case, hence Lemma 1-5 apply. For example, the following 
Hamiltonian (with qpp = (0, 0, 1, 0)) satisfies the genericity assumptions (for 
MO): 

Hoix,v,D;c) = {^ + b)^{D-lY-cD + ^v^ + ^arx^l-D) + ^a2X^ + ^a3X* 

(3.3.3) 

Taking as > 0, 6 > corresponds to compact energy surfaces Ho{x, v, D; c) 

h for any finite h. Adding the perturbation Hi = e(l — ^x"^) cos{k9), the cor- 
responding vector field is of the form: 

X = V (3.3.4) 
V = ai{D — l)x -\- a2X^ -\- a^x^ -\- ex cos{k9) (3.3.5) 

e = + b){D - 1) - c - l-aix^ (3.3.6) 

2(23 2 

i) = ke{l - ^x^)sm{k9), (3.3.7) 
attaining hyperbolic resonances for c > at: 

Dr{c) = 1 + (3.3.8) 

2a3 + 

a parabolic resonance at Dp = 1, Cp = and a fiat parabolic resonance 
at these values for 6 = 0. Numerical experiments with this model produce 
similar results to the ones presented for the atmospherical problem described 
in section 2. 
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4 Conclusions 



The behavior near parabohc resonances provides a new type of chaotic mech- 
anism for two d.o.f. Hamiltonian system. It has been estabhshed that this 
phenomenon is of co-dimension 1, hence it is expected to appear in numerous 
apphcations. Moreover, the basic idea of Lerman and Umanskii [15], that the 
addition of another non-separable d.o.f. may be considered from bifurcation 
theory point of view as the addition of another parameter, applies to higher 
degrees of freedom as well - e.g. we may think of the wave speed c as another 
conserved quantity of a larger system showing that for higher dimensional 
systems the parabolic resonance case may be generic. 

Preliminary numerical observations of the behavior near parabolic and 
flat parabolic resonances are presented in figures 2.0.1 - 2.0.5. The features 
of the parabolic resonance seems to mix three different dynamical behaviors - 
localization as in the elliptic resonance case, homoclinic chaos which destroys 
this localization and slow passage through a bifurcation which leads to strong 
instabilities. 

For the atmospherical model with {B[(f>) = 0), for any given wave speed 
— 0.5<c<0,a hyperbolic resonance occurs: particles sent from the equa- 
tor's vicinity with westward velocity Uq c(< 0) participate in a complex 
dynamical behavior moving north and south chaotically, reaching latitudes of 
arccos(l — 2|mo|) + 0{^/€). The motion of such particles may be non-uniform 
zonally [1] (i.e. some zones will be visited more frequently than others), yet 
generally, zonally localized motion is not robust. When the speed of the 
travelling pressure wave vanishes (i.e. the perturbation is of the form of a 
standing wave) parabolic circles of fixed points occurs for e = 0. Small per- 
turbations of the form of standing waves or travelling waves with small wave 
speeds results in parabolic resonance. Moreover, this case is degenerate lead- 
ing to strong instabilities as shown in figure 2.0.1. This degenerate model 
thus supplies a mechanism for the transport of particles with small initial 
velocities near the equator to high latitudes. Adding a latitude dependent 
pressure gradient {B[(f>) ^ 0), which physically corresponds to including the 
infiuence of jets [1] removes the degeneracy. If this gradient is not large 
(which seems to be a reasonable assumption from the physical point of view 
as the pressure gradient is proportional to the particles velocity, which should 
be of order 0.01) the degenerate instabilities still appear (figure 2.0.3). 

The results presented are preliminary - analytical and extensive numerical 
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studies are under progress. It appears that near a generic parabolic resonance 
the instabihty zone in the latitude coordinate corresponds approximately 
to the maximal latitude achieved by an unperturbed homoclinic orbit to a 
hyperbolic circle with D = Dp — 0{^/€). The behavior near a flat parabolic 
resonance is not resolved yet even numerically. In particular, the scaling of 
the radius of instability with respect to the perturbation strength (if it is 
finite), and with respect to the distance in phase-space and parameter space 
from the fiat parabolic resonance case are yet to be found. General conditions 
for the occurrence of a fiat parabolic resonance are yet to be derived. 

Many questions are left unanswered at this point, and in particular the 
relations of this work to McLaughlin, Klafter, Meiss, Melnikov and Zaslavsky 
presentations are yet to be explored. 
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